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Turbulence in the interstellar medium has been an active field of research 
in the last decade. Numerical simulations are the tool of choice in most cases. 
But while there are a number of simulations on the market some questions 
have not been answered finally. In this paper we are going to examine the 
influence of compressible and incompressible driving on the evolution of tur- 
bulent spectra in a number of possible interstellar medium scenarios. We 
conclude that the driving not only has an influence on the ratio of compress- 
ible to incompressible component but also on the anisotropy of turbulence. 



1. Introduction 



Turbulence in the context of the interstellar medium (ISM) was for the first time dis- 



cussed by von Weizsacker (1951). This discussion, however, was soon abandoned by the 



scientific community. Until the late seventies there was not much progress in the field of 
interstellar turbulence. Only when the observational techniques became better at this 
time it became clear that there were power law correlations between different structures 
in the ISM. At first turbulence was only connected to the smallest scales of the ISM but 
with improved observations it was found that turbulence was present at all spatial scales. 



A corresponding spectrum of the density fluctuations can, e.g., be found in Armstrong 



et al. (1995) 



Besides the well known density fluctuation spectra also velocity power spectra may 
be obtained. One possible method for this has been presented by |Lazarian fe Pogosyan| 
(1997) and was applied by Padoan et al. (2006). The connection between the velocity 
and density observations are not fully understood, as pointed out by Klessen (2000). For 



the case of the magnetic field structure functions may be observed with radio telescopes 



(e.g., Haverkorn et al., 2004). For a more detailed account of the observational evidence 



for interstellar turbulence see Elmegreen Sz Scalo (2004). 
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Since there are yet no in-situ measurements of the interstellar medium, a number of 
questions is still not resolved. One of those is the question of compressibility. Many 
theories assume an incompressible medium, yet this is not compatible with the density 



fluctuation observations in the ISM. In the wave picture of interstellar turbulence (Lazar 



et al. , 2003 Spanier & Schlickeiser 2005 ) it may well be assumed that up to certain wave 



numbers a major compressible component is present. Another issue which is not yet fully 
answered is the question of the driving of the turbulence. The energy balance of the 
interstellar medium suggests that supernovae are the driving force and energy is dissi- 
pated at small scales. Whether this driving is compressible or incompressible is not clear. 
On the one hand the driving scenario involving supernova shock waves is connected to 
compression, since the shocks themselves are compressed. On the other hand streaming 
cosmic rays may be also a source of driving especially through incompressible Alfven 
waves (Skilling, 1975). 

In the course of this paper we will come back to the main theories for anisotropic MHD 
turbulence and review recent simulations. We then present the results of our simulations 
and their implications for the understanding of compressible MHD turbulence. 

In this paper we will at first discuss the current knowledge about the theory of plasma 
turbulence, then we will p resent our numerical methods. The results are shown for 
short time scales of turbulent driving, the evolving spectru m during driving and the 
saturated spectrum. 



2. Theory 



2.1. Turbulence theory 

From the theoretical side the discussion about a turbulence theory started out from 
the investigation of incompressible turbulence described by the Navier-Strokes equation. 
The early discussions by Kolmogorov (see, e.g., Kolmogorov, 1941) used the assumption 



of spatially homogeneous turbulence. Later on, however, it was found that intermittency 
had to be taken into account leading to a slightly different form of the velocity spectrum. 



For different intermittency models and a general overview of the turbulence see Frisch 



( 1995 ). For the turbulence in a magnetised fluid the model by She & Leveque ( 1994) was 



extended to include compressibility and a magnetic field. Politano & Pouquet (1995) 



have discussed the role intermittency in the solar wind, proposing sheet like structures 



(also discussed in Grauer et al. 1994). 



The next important step is the inclusion of a background magnetic field. In Shebalin et al. 



(1983) numerical studies have been conducted showing anisotropics of MHD plasmas 



due to the background field. This aspect was addressed by Sridhar h Goldreich (1994), 



where the authors use a heuristic wave-model to get some idea of the anisotropy of 
the turbulence spectrum of magnetic turbulence with a background field. The effects 
of a background magnetic field, however, are by no means fully understood. This is 



reflected by the debate started in Goldreich & Sridhar (1995), where a discussion about 



the correct description of parallel and perpendicular cascades started. Another aspect is 
the question when turbulence can be described by a wave-interaction picture and when 
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non-linear processes dominated instead. 

To resolve this question a number of simulations have been conducted in the last 
15 years. The first numerical test has been conducted by Cho & Vishniac (2000), a 



prominent example are the simulations by Maron & Goldreich (2001), who consider 



simulations of decaying incompressible turbulence. The incompressible plasma attracts 
far more attention by modellers for a number of reasons: A prominent one is that 
incompressible plasmas may be simulated with spectral methods providing a spectrum 
extending to high wave numbers. Apart from that basic turbulence theories are mainly 



developed for incompressible fluids. But also Chandran (2005) has done additional work 
on compressible turbulence in the ^ < 1 limit for the evolved spectra. Due to the 
very low (3 limit sound waves are neglected in these simulations. Here, however, we will 
show the importance of sound waves in the context of compressible plasmas. A work 



which is more closely related to ours is the one of Kuznetsov (2001), where the authors 



investigate compressible plasmas with high and low values of /3. In his approach only 
weak turbulence is considered. As mentioned above it is not yet clear if a wave-picture 
is applicable, but in the weak turbulence limit the wave-picture is explicitly assumed. 
Additionally a theoretical description including major cooling and damping mechanisms 
for the interstellar medium has been presented by Lithwick & Goldreich (2001). Here 
different (5 regimes are discussed. 

2.2. Excitation of MHD waves with compressible or incompressible 
distortions 



In Federrath et al. (2008, 2010) a discussion similar to our discussion regarding the influ- 



ence of compressible versus incompressible driving has been started. But this discussion 
was limited to hydrodynamic simulations of the interstellar medium. A more recent pa- 



per (Federrath et al. , 2011 ) added the magnetic field and discussed the effect of turbulent 



magnetic field growth. We will extend this discussion in the case of magnetised fluids to 
the anisotropics present in the driven turbulence. This discussion has been continued by 



Lemaster & Stone (2009), who did simulations of MHD plasmas, but limited themselves 



to incompressible driving but with variations of the peak wave number of the driver. In 
their work it was shown that approximately a factor of 128 between the driving peak 
and the total resolution is sufficient for convergence. 



2.2.1. Compressible driving 

We inject compressible disturbances isotropically into our simulation domain. These 
disturbances can distort the magnetic field and so excite waves. The induction equation 
shows us which waves can be expected. Without loss of generality the preferred direction 
of the magnetic field can be chosen to be the x-direction (defined by the unit vector e x ). 
For the induction equation we find that: 



dB 
~dt 



-V x E = V x (Sv x Be x ) 



3 



Where we used E = —5v x B for the electric field which holds for ideal plasmas. If 
the compressible fluctuations are in x-direction no magnetic field will be injected as 
(5ve x x B ■ e x ) = 0. If the distortion is perpendicular to the background field (e.g. 
5v = ve y ) this results in a distortion of the magnetic field in x-direction: 



ike y x (5ve y x Be x ) (2) 



dB 

Ik 

dB . 

— = ike y x (-v)Be z 

OB 

— = -ie x kvB (3) 

Where we used V = ike y as for compressible waves k is parallel to 5v 

For the case Sv = ve z we also find a distortion of the magnetic field in x-direction: 

dB 

-— = ie x kvB (4) 
dt w 

As the distortion of the magnetic field is perpendicular to k, we would expect to find 
compressible waves that propagate in perpendicular direction with respect to the mean 
magnetic field. We can also identify these waves as for MHD-waves the propagation in 
the perpendicular direction is only possible for fast magnetosonic waves. 

2.2.2. Incompressible driving 

Here we want to figure out if Alfven waves being the only incompressible waves in MHD- 
plasmas can be generated form those incompressible fluctuations. As for incompressible 
fluctuations k is parallel to the background field V = ke x . So we find for the magnetic 
field: 

dB 

-7^- = ke x x (5ve y x Be x ) = e y kvB (5) 

or 

dB 

— = ke x x (5ve z x Be x ) = e z kvB (6) 

As the distortion of the magnetic field is perpendicular to k this results in Alfven waves 
propagating parallel with respect to the mean magnetic field. 

2.3. Three wave interaction 

Three wave interaction is used to describe the interaction of compressible and incom- 
pressible waves in a weakly turbulent plasma and has been described in the first place 



by 


Chin k Wentzel 


(1972 


Vainio & Spanier 




!005 


)• 



For our investigations where we either drive with pure compressible or pure incompress- 
ible fluctuations one has to put the question with which process compressible fluctuations 
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are generated from incompressible fluctuations and the other way around. As this hap- 
pens for small turbulent amplitudes we claim that three wave interactions are the key 
process. 

Here a validation is carried out how Alfven waves are generated by fast magnetosonic 
waves. It is also shown what can be predicted about the propagation direction of the 
resulting Alfven waves. The calculation of the generation of fast magnetosonic waves 
from Alfven waves can be done analog. 

In principle Alfven waves can be generated from fast magnetosonic waves by one of these 
two processes: 

A + + F- <- A ± (7) 
A+ + A- -> F ± (8) 

Here we will restrict the discussion to the special case of perpendicular propagating fast 
magnetosonic waves with (3 <C 1. For this we find the phase velocity of these waves to 
be: 

v ph « {y/T+P)v A (9) 

In this case the reaction equation Eq. [9] provides no valid solution when computing the 
resonance condition (i.e., the energy equation). For reaction equation [8] , however, we 
find the following resonance condition: 

{ ) = v A (y/l + P)k$ 

v A (y/T+P)((k+-k2..) 2 + l4 



*M(fe£|| + k A,\ 



v a{K h + k A ,,) 



a/2 



(10) 
(11) 



with 



k: 



+ K 



(12) 
and 

so with equation 11 and 12 we find for the parallel component of the fast magnetosonic 
waves: 



n A,\\ All ~~ -^ n F,\\ n A,± ^ A,. 
As the fast magnetosonic waves propagate in the perpendicular direction kp\ 



k 



A,\\ 



2**11 



For the perpendicular component of the fast magnetosonic waves we find 

1 



-kp 



With that we find that 



tan ( 



(13) 



(14) 



(15) 



k Al x \k F 

where 6 is the angle between the propagation direction of the fast magnetosonic waves 
and the Alfven waves. As /? <C 1 it can be followed that tan0 ~ 1. This gives us a 
preferred propagation direction for the Alfven waves close to 9 = 45 degree. 



A similar study has been undertaken by Chandran (2008) for the low-/? case taking 



into account also fast magnetosonic wave. Lithwick & Goldreich (2001) have discussed 



systems with high and low j3, their approach is an analytical one. 
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3. Numerical Simulations 



3.1. MHD equations 

For our studies we used the ideal MHD equations in a periodic domain, where the system 
is closed by an isothermal equation of state. Since our analysis is a principle one we will 
give all variables in non-dimensional, normalised units. For that purpose all variables 
are split into a dimensionless variable of order unity and a normalisation constant. Here 
we use four independent variables for the normalisation: The length of the simulation 
domain L, the mass of the hydrogen atom mo, a typical number density no, and the 
temperature of the system To which directly relates to the speed of sound c s . Therefore 
length scales are normalised as x = Lx, velocity as u = uqu = c s u, density as p = mono/5 
and so forth. The resulting set of normalised MHD-equations is the following: 



dp 
dt 



(16) 



£ = _v(f + (p + f)i-ai) + * < 17) 

p = ~p = p 19 
Here p is the mass density, s = pvis the momentum density and B indicates the magnetic 



induction and p is the thermal pressure. The isothermal equation of state - Eq. ( 19 ) - 
is simplified due to our choice for the normalisation constants. 

In the following we will classify our results on the basis of f3 we used for our simulations. 
We only use it for the initialisation of the unperturbed background magnetic fields. It 
is defined as the fraction of the thermal and magnetic pressure. 

3.2. Numerical model 

For the time evolution of the MHD-equations we use a second order CWENO (Centrally 



Weighted Essentially Non-Oscillatory) algorithm (see e.g. (Kurganov &; Levy, 2000)). It 
is very convenient for the simulation of shocks as we find them in fluids and plasmas. In 
the mathematical description these shocks are so called "Riemann Problems" and their 
accurate numerical treatment with a Godunov solver results in high numerical costs. 
The CWENO scheme avoids to solve these Riemann Problems directly by averaging 
over the appropriate fractions of every cell. Additionally oscillations are suppressed by 
the reconstruction. This results in a fast scheme for the treatment of the MHD equations 
with passable small dissipation. This CWENO method is combined with a third order 
Runge-Kutta algorithm which has been chosen because of the small memory costs with 
comparatively high resolution in time, due to the fact that only one additional field has 
to be buffered. 
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Besides the direct numerical simulation (DNS) technique also the large eddy simula- 
tion. This is for example demonstrated in Chernyshov et al. (2007, 2009) for the case 



of decaying turbulence. We are refraining from using this method since it requires a 
number of implications on the turbulent cascade on small scales which we are not sure 
of. 



3.3. Turbulence driver 

For the external stirring of the fluctuations we define a function / in Fourier space with 

f k = s • AT 7 / 4 exp(2vrzp) (20) 

where "s" and "p" are random numbers between zero and one. The random number 
generator for the former of these obeys a Gaussian distribution, whereas the latter is 
uniformly distributed. This function yields the Fourier space distribution of the 
velocity fields used for driving the fluctuations in our simulations. After transforming to 
configuration space fk — > f x we make the distinction between the cases of compressible 
and incompressible driving. For the former we use 5v = V/ x , whereas for the latter we 
employ 5v = V x / . As / has, thus, to be a vectorial function in the case of incompressible 
driving we define each component separately with different random numbers. Each 
component of the wave- vector k in our Fourier space spectra ranges from one to eight 
yielding a maximum absolute value of | k |< 13. The perturbations are set up in a form 
such that no net momentum is added to the numerical domain. The energy input is 
normalised to an estimate for the average energy input by supernovae into the ISM of 
our Galaxy into a volume of the size of the simulation domain. The resulting driving 
spectrum is isotropic and it yields v 2 (k) oc A; -3 / 2 which yields a production spectrum 
harder than the anticipated inertial spectrum. 

Finally we define the compressible and incompressible energy spectra. Both are cal- 
culated in Fourier space on the basis of the velocity fields with 



comp 



P 



mcorap 



(k) 
(k) 



1 



k ■ v(k) 
k x v(k) 



(21) 
(22) 



where k is the normalised wave vector. The overall energy spectrum of the velocity fields 
is 



1 



PK(k) = -v\k) 



(23) 



where the sum of P C omp{k) and Pincomp(k) is Px(k). For the following discussion we have 
to make an important differentiation: Pk is the power spectrum of velocity fields in the 
whole system, whereas P^ is the power spectrum of the velocity fields of the driver. 
This is important because the turbulent energy of the velocity fields of the whole system 



7 



doesn't have to accord to the driving spectrum. This will be important in chapter 4.1.2 
Additionally we define the omnidirectional energy spectrum of the magnetic fields with 



Pb = Y p tB 2 (k). (24) 



3.4. Simulation setup 



For all our investigations we use in principle the same simulation setup: We start out 
from an unperturbed plasma with a homogeneous background magnetic field and inject 
with each time step either compressible or incompressible energy with our turbulence 



driver explained in section 3.3 We do this for compressible and also for incompressible 
driving and both runs are performed for /3=0.01 and /3=10. This enables us to analyse 
the influence of the different drivings on as well the highly magnetised as the quite 
hydrodynamic plasmas. So the basis of the following research are those four runs. With 
the simulation output we gain during the simulation we calculate the turbulent energy 



spectra defined in section 3.3 



4. Results 



Here we will discuss three distinct topics relating to our simulation results: First we 
investigate the physics on short and intermediate timescales where the fluctuations of 
the velocity fields excite modes. This is done in section 4.1. Then the time evolution of 
the energies is studied in section 4.2 Finally an investigation of the convergence range 
is done in section 4.3 In order to quantify the anisotropy we use a scalar parameter 



which is derived from 2D spectra. The parameter is described in detail in [A} 



4.1. Short time scales 

This investigation is of special interest because it can be observed how the first waves 
are generated for compressible as well as for incompressible driving and how these waves 
interact with each other. On the basis of this research we can determine which pro- 
cesses dominate the turbulence for small amplitudes of the fluctuations. When taking 
the supernova injection scenario for granted, the build-up of turbulence is a transient 
phenomenon which may even be observed close in time and space to the energy source. 
One observable could be accelerated particles in the upstream of supernova and their 
respective electromagnetic emission. 

For the investigation of the two dimensional spectra we need to keep in mind that 
after this short amount of time we have to distinguish between the driving region that 
is dominated by the driving spectrum injected for each time step and the part of the 
spectrum that can only be reached via a cascade. After this short amount of time the 
cascading process has only taken place to some extent so we concentrate mostly on the 
driving region. Via the intersection points of the contour lines with the k±- and fey -axes 
we can access the anisotropy of the turbulent fields. 
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4.1.1. P B 



Figures 1(a) to 1(d) show the two dimensional spectra of Pg. We see for both drivings 
and values of the plasma-/? a very anisotropic spectrum (/c^ mso (l) = 2) within the driving 
range. The spectra show a remarkable feature: For incompressible driving contour lines 
in the driving range rise from the k± axis (following the diagonal) , while for compressible 
driving the contour lines rise from the fcn axis. The slope of the contour lines indicates 
a spectrum which has in parts a positive exponent with respect to k. 

There is an explanation for this behaviour, which may be deduced from the fact, 
that the difference exist between compressible and incompressible fluctuations. For the 
case of incompressible fluctuations, the vector 5v is perpendicular to the wave vector k, 
while for compressible fluctuations, the both vectors are aligned. If we examine now the 
induction equation for incompressible fluctuations (Eq. [5J moving in parallel direction, 
we find 

— = V x (5ve y x B e x ) (25) 
this yields for the turbulent velocity spectrum 



(26) 



This yields a rising spectrum in parallel spectrum as seen above. It also explains the 
crossing of the contour line and the k± axis. For the compressible driving the explanation 
follows the same line, but is essentially different, since the vector 5v has a different 
orientation. 

We can make further distinction between f3 = 0.01 and f3 = 10: If we examine the 
2D spectra closely we find a cutoff in some of the spectra, which may be identified by 



very close contour lines (as for example in Fig. 1(a)). For compressible driving we find 
a strong cutoff for /3 = 10, while it seems to have dissolved already for /3 = 0.01. Since 
the induction equation makes no distinction between the two values of f5 it seems that 
another physical process is taking place here. By examining Pcomp and Pshear we will 
show later that this may be caused by a mixture of Alfven and fast magnetosonic waves. 
2(a) to 2(d) the magnetic field after a normalised time of 5 • 10 -3 can be seen. One 



In 



may now compare the eddy size and spectrum (as done for example in Beresnyak Sz 



Lazarian|2009| ) to see at least qualitatively that the elongation of the eddies is correlated 
with the cutoff seen in the two-dimensional spectra. This is of course not unexpected 



4.1.2. P K 

Figure |3(a)| to [3(d)| show the two dimensional spectra of Pk for /3=0.01 and /3=10 for 
both compressible and incompressible driving. 

For /3=0.01 and incompressible driving we see in the driving region of Pk a perpendicular 
preferred direction (/c^ mso (1.5) = 1.8) whereas for /3=0.01 and compressible driving we 
find a parallel preferred direction (Ar^ mso (1.5) = 1.5) although P K is for both cases 
isotropic. The term "preferred direction" is used here when the two dimensional spectra 



9 




in /cm In/cy 
(a) /3=0.01 incompressible driving (b) /3=0.01 compressible driving 




m/cii In/cy 
(c) ,3=10 incompressible driving (d) 13=10 compressible driving 



Figure 1: Turbulence spectra Pb after t = 5 • 10 3 simulated on a 512 3 grid. The color 
denote the natural logarithm of the numerical energy. 
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X X 
(c) /3=10 incompressible driving (d) /3=10 compressible driving 

Figure 2: Turbulent magnetic fields after t = 5 • 10~ 3 in real space on a 512 3 grid. 
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(a) /3=0.01 incompressible driving 
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(c) /3=10 incompressible driving 
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in Am 



(b) ,3=0.01 compressible driving 




In A" 1 1 

(d) /3=10 compressible driving 



Figure 3: Turbulence spectra Pk after t = 5 ■ 10 3 simulated on a 512 3 grid. 



are anisotropic. The contour lines of a spectrum with a parallel preferred direction 
cross the axes for higher values of fey than k±. For the contour lines of a spectrum 
with a perpendicular preferred direction it's the other way around. For /3=10 we see for 
compressible as well as for incompressible driving a very homogeneous driving spectrum. 
This can be understood quite easily. We drive the turbulence with compressible or 
incompressible fluctuations in the velocity fields. As shown in chapter 4.1.1 this results 



in an excitation of MHD- waves in the driving range, so turbulent energy of the velocity 
fields is transferred into magnetic energy. In the direction of preferred excitation of 
waves this energy is missing in the velocity fields and so in Pk- This is why for the case 
of incompressible driving and (3 = 0.01 we find Alfven waves propagating parallel to the 
mean magnetic field and so a perpendicular preferred direction in Pk- For compressible 
driving we excite fast magnetosonic waves perpendicular to the mean magnetic field 
which results in a parallel preferred direction of Pk- 

For (3 = 10 although the same is in all likelihood happening this effect can't be seen, 
because the magnetic fields are much weaker and so less waves are excited and through 
this less energy of Pk is transformed into Pg. Hereby the missing energy in Pk in one 
of the preferred directions can't be detected in Pk- 
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4.1.3. Pshear and P C 



omp 



Now we will investigate the two dimensional spectra of Pcomp and Pshear- Figure 4(a 
to 5(d) show the two dimensional spectra of Pcomp and Pshear at a normalised time of 
5 • 10" a . At this very early time the type of driving dominates the turbulent energy 
of the velocity fields. So for j3 = 0.01 the driving component shows of cause the same 
behaviour as Pk for f3 = 0.01. We also find for compressible driving a parallel pre- 
ferred direction in Pcomp (A;^ niso (2.2) = 2.0) as in this case fast magnetosonic waves are 
generated propagation perpendicular to the mean magnetic field and this energy comes 
from the turbulent energy of the velocity components. For Pshear and incompressible 
driving we find a perpendicular preferred direction {k^ lso {\A) = 1.8), as we have an 
excitation of parallel propagating Alfven waves. The non driven components - Pcomp for 
incompressible driving and Pshear f° r compressible driving - show in principle the same 



diagonal symmetry. In chapter 2.3 it was shown that Alven waves can be generated via 
three wave interaction from fast magnetosonic waves propagating perpendicular to the 
mean magnetic field and that these Alven waves propagate with an angel of 45 degree 
to the fast magnetosonic waves. This prediction corresponds to the two dimensional 
spectrum of Pshear for compressible driving. The same can be done for Pcomp and in- 
compressible driving. 

For (5 = 10 we find the same Pcomp and Pshear limitations as for Pk'. Because of the 
weak magnetic fields only few waves are generated and this effect becomes as little no- 
ticeable as in the driven component as in Pk- As the non-driven component is also 
absolutely isotropic, no clear conclusion about the dominating physical process can be 
made. 

One should note that there should exist a connection between the P CO mp results and 
the density fluctuations through the continuity equation. Density fluctuations have been 



studied in detail by Kowal et al. (2007), where only solenoidal driving is used. 



4.2. Evolving spectrum 

In the following chapter Pr ,Pc amp ,P shear and Pb for the four different setups are plotted 
against the normalised time (see figure 6(a) and |6(b) as well as figure [8(a) and [8(b)] 



It can be seen that there is no turbulent energy in the fields at the beginning of the 
simulation. As the turbulence is driven continuously through the whole simulation more 
and more turbulent kinetic energy is injected in the fields with time until at last the 
turbulence saturates (this is when the input energy from the driving is balanced by the 
numerical dissipation) . A detailed analysis about the saturation range is given in section 



4.3 Here we are mainly concerned about the physics before convergence is reached. The 
absolute amplitude of the energies in the convergence range is influenced by numerical 
effects and is so neglected in the discussion. 

4.2.1. Turbulent kinetic energy Pk 

The temporal evolution of Pk differs fundamentally between compressible and incom- 
pressible driving (see figure [6(b) ) . This is especially apparent for f3 = 10. For com- 
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(a) /?=0.01 incompressible driving 
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(b) /3—Q.01 compressible driving 
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(c) (8=10 incompressible driving (d) /3=10 compressible driving 

Figure 5: Pshear after t = 5 • 10~ 3 simulated on a 512 3 grid. 
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Figure 6: Turbulence spectra Pk and Pb as a function of the simulation time 
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pressible driving we find a distinct maximum for Pk at t~ 0.05. After that the energy 
decreases before finally saturation is reached. 

What causes the formation of this maximum? In principle we find for both - compress- 
ible and incompressible fluctuations - the same behaviour: If the magnetic field lines are 
bent, the cascade and so the dissipation of the energy in the system is more efficient. 
From the data we can see that for compressible driving perpendicular fluctuations are 
generated that decay into the parallel direction which leads to an isotropic spectrum. 
To do this they need a component parallel to the background magnetic field. If the 
background field is quite tangled the fluctuations have a parallel component locally and 
can decay. For incompressible driving parallel fluctuations are generated that decay into 
perpendicular direction and so they need a perpendicular component to the background 
field. This is why they as well can only decay in a tangled field. The formation of the 
maximum can occur if the driving energy will bend the magnetic field lines only insuffi- 
ciently. So for quite a long time there won't be a cascade and the resulting dissipation 
of the energy. Hereby lots of energy is accumulated until finally the field lines are bent 
enough to build up the cascade. Then the accumulated energy for small k can at last 
be reduced by the cascade and the combined dissipation. If the driving energy bends 
the magnetic field lines efficiently from the beginning the cascade will build up from the 
beginning. Through this no energy is accumulated and no maximum will occur. 
The magnetic field lines for compressible and incompressible driving at this maximum 
(t=0.05) and after convergence is reached are shown in figure |10(a) to 10(d) It can 



be seen that for the case of compressible driving the magnetic field lines are still quite 
smooth at the maximum whereas in the convergence range they are quite tangled. For 
the case of incompressible driving however the field lines are already quite tangled at 
t=0.05. This can be understood, when considering how the driving acts on the fluid. 
For a incompressible driving we are basically adding eddy-like motions to our fluid. This 
means that we get some shear-flows, which will certainly bend the background field lines. 
When adding a compressible velocity field, however, the background field even remains 
undisturbed for some of the input wave-modes. This shows, that an incompressible driv- 
ing field will tangle the field lines much more efficiently Thus, we find for (3 = 10 that 
with the incompressible driving field the field is tangled so early that such an energy 
storage as in the case of the compressible driving can not form at all. 

This tangling of the magnetic field lines can also be accessed quantitatively. For this 
we computed the spatial averages of the quantity 

a = 1 - 2/u 2 (27) 

where \x is the angle between the local and the initial magnetic field. For this Parameter 
we find alignment along the initial magnetic field direction for a = — 1, where a magnetic 
field that is perpendicular to the initial one would yield a = 1. Due to the fact that there 
is a higher degree of freedom for the perpendicular direction an isotropic distribution 
would yield a = 1/3. With this in mind we can discuss the results for our simulations. 

For an initial plasma-/? of j3 = 0.01 the alignment-parameter a even for the saturated 
turbulence at t = 0.3 never exceeds values of -0.9, which is still very parallel. For (3 = 10, 
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Figure 7: Distribution function of the magnetic field angle for different times and differ- 
ent numerical simulation runs - all results are shown for /3 = 10. On the left 
we show results for the incompressibly driven turbulence for t = 0.01 (black 
curve) and for t = 0.3 (red curve). On the right the results for compressible 
(black curve) and incompressible driving (red curve) are compared for t = 0.05. 
In both plots the dashed line indicates homogeneous distribution. 



however, the situation is very different. In this case we find at t = 0.05 (which is the 
time, when the peak in the energy evolution is most prominent) on the one hand a nearly 
isotropic distribution with a = 0.296 for the incompressible driver. For compressible 
driving on the other hand we find a = —0.158576, where both of these simulations show 
isotropic distributions at t = 0.3. This shows that the compressibly driven simulations 
with j3 = 10 are the only ones, for which the tangling of the fieldlines still changes 
significantly between t = 0.05 and the saturated state. This is also illustrated by the 
distribution functions shown in Fig. [7] for some of the simulations with /3 = 10. On the 
left we show the distribution function for very early times, where the magnetic field is 
still pretty much aligned, and the distribution for the saturated state, which is obviously 
quite isotropic (here one has to take into account that an isotropic distribution in all 
spatial directions yields a distribution of the form sin# when projected from 3D to a 
Si I — B± distribution- function) . On the right we compare the situation for incompressible 
and compressible driving at the time of the maximum. Apparently the distribution is 
still much more aligned for the compressibly driven simulation. 

In this context the simulations of Beresnyak &; Lazarian ( 2009 ) should be mentioned: 
Taking into account a partially similar situation, they evaluate also multidimensional 
spectra, but they follow a different approach. While in this paper the global frame of 
reference is used (the axes are determined by the background magnetic field), they are 
using a local frame of reference. The discussion which of these two approaches is the best 
is tedious and maybe not helpful. Looking at Fig. [7] we see that there is a clear bending 
of field lines, which would suggest using the local frame. But in our opinion there are 



17 



two major drawbacks using the local frame: i) the definition of the local frame is not 
completely unique and may be governed by cutoff effects, ii) in the case of very weak 
turbulence local and global frame will be equivalent, while for very strong turbulence we 
may see differences in the local frame picture, which will in the end not affect the global 
picture (which is usually the one being observed). 



4.2.2. P B 

For early times it can be seen that Pb has higher values for compressible than for 
incompressible driving (see figure 6(a)| and table [T]) . This is due to the fact that for in- 



compressible driving Alfven waves are generated propagating in the x-direction whereas 
for compressible driving fast magnetosonic waves are generated propagating within the 
whole perpendicular y-z-plane. So the magnetosonic waves have an additional degree 
of freedom for their excitation and because of that more fast magnetosonic waves than 
Alfven waves are excited for the same driving strength and driving time. 
During the whole simulation much more turbulent energy of the velocity fields is con- 
verted into magnetic energy for small (3 than for high f3. This is due to the fact that the 
background magnetic field is stronger for small plasma- /3s which yields stronger mag- 
netic fluctuations by tangling the velocity fields. 

Pb sometimes also shows a maximum before it reaches the saturation range (e.g. for 



(3 = 0.01 and incompressible driving at 0.09 in figure 6(a)). But in this case it barely 
depends on the driving but on (3. This is because the rate magnetic energy is converted 
from the velocity fields depends strongly on (3. For (3 = 10 this rate is anyway so small 
that even a quite inefficient cascade is adequate to remove enough energy that no max- 
imum can build up. For small f3 the energy is converted fast enough that in principle a 
maximum can build up. As for the case of incompressible driving more energy is con- 
verted in the saturation range the maximum is even slightly higher (see figure 6(a)| for 



(3 = 0.01 and incompressible driving at 0.09 as well as for /3 = 0.01 and compressible 
driving at 0.075 ). 



4.2.3. Pshear and P c 



ornp 



Pshear and Pcomp are plotted in 8 (a 



9(d) together with the total 



to 8(b)| and |9(a)| to 

energy of the velocity fields Pk- In table 2] the energies are listed for a normalised time 
of 5 • 10~ 3 and for the saturated state of the spectrum. 

At the beginning of the simulation the energy of the driving component is higher for 
(3 = 10 than for f3 = 0.01. This is because for bigger plasma- /3s there are smaller 
magnetic fields and consequently less waves are generated. So less energy of the driving 
component is transformed into magnetic energy and hence missing in the energy of the 
driving component. Additionally as there are less excited modes for /3 = 10 the non 
driven component is smaller because less modes are generated that can decay into other 
waves and so transform compressible modes into incompressible modes and the other 
way around. 

Whereas the driving dominates Pk in the beginning it stands out that the incompressible 
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Figure 8: Turbulence spectra Pcomp and Pshear as a function of the simulation time 
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Table 1: magnetisation of the system. 



energy after some time dominates Pk independently of the kind of driving and the 
plasma- (3. The dominance of the incompressible energy might be caused by a more 
efficient conversion of compressible modes into incompressible modes than the other way 
around. The saturated results given in chapter |4.3| support this idea. 



4.3. Saturated state of the spectrum 

The saturated state of the spectrum has been described in a number of previous papers. 



In the work of Cho & Lazarian (2002, 2003). While Cho & Lazarian (2002) concentrated 



on the simulation of low-/? plasmas, Cho & Lazarian (2003) also considered high-/3 plas 



mas. In both cases incompressible driving was used. As discussed above we will also 
use compressible driving, which may be motivated by supernova shock waves driving the 
ISM. Another point which is different is that the authors of the aforementioned studies 
attempted the decomposition into different wave modes. Our approach will focus on the 
decomposition in irrotational and solenoidal fields. This is specifically useful when one 
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Figure 9: Total turbulent energies Px^Pcomp and Pshear as a function of the simulation 
time 
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Table 2: compressible and incompressible energy of the system. 
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(a) magnetic field lines for /3 = 10 with (b) magnetic field lines for j3 = 10 
compressible driving after t=0.05 with compressible driving after t=0.3 




(c) magnetic field lines for /3 = 10 with (d) magnetic field lines for f3 — 
incompressible driving after t=0.05 10 with incompressible driving after 

t=0.3 



Figure 10: Magnetic field lines for /3 = 10 with compressible and incompressible driving 
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assumes that the turbulence is dominated by shocks. 

Here we present the two dimensional spectra for the saturated state of the turbulence. 
For f3 = 0.01 and compressible driving Pk, Pshear and Pcomp have a clearly parallel 
preferred direction (see figure 12(b)|13(b) and |14(b)[ ). In the saturated state we are 
concerned with fully nonlinear interactions, which form an anisotropic cascade, as we 
see from the results. For incompressible driving we find that Pk and Pshear develop a 
perpendicular preferred direction (see figure 12(a)| and 13(a)| ) , as it has been predicted 
by Shebalin et al. (1983). The spectrum of Pcomp however is quite isotropic (see figure 



14(a)). As it has been already mentioned in chapter 4.2.3 we believe that compress- 



ible modes may easily be converted into incompressible modes, but not the other way 
around. This is confirmed by the two dimensional spectra, because for compressible 
driving P shear has the same anisotropy as Pcomp- Actually the incompressible modes 
would build up a perpendicular cascade as can be seen by the results in 12(a)| and 13(a 



as they are however directly generated by the compressible modes, the cascade is forced 
to build up the same preferred direction. Contrariwise in the case of incompressible 
driving Pcomp has not the same preferred direction as Pshear- So the transformation 
of compressible modes from incompressible modes is not that efficient. This is why a 
superposition of compressible and incompressible modes and so of a parallel and per- 
pendicular cascade can be seen. The result is an isotropic spectrum. 
For (5 = 0.01 it can be seen that for incompressible driving Pg has the same perpendic- 
ular preferred direction that we also find for Pjc, Pshear (see f igure |ll(a) ) whereas for 
compressible driving there is no preferred direction (see figure |ll(b)| ). In figure |9(a)| it 
can be seen that for incompressible driving much more incompressible modes than com- 
pressible modes are in the system. This is why we would also expect a perpendicular 
cascade to build up. For the case of compressible driving (see figure [9(b)] ) we have only 
slightly more incompressible than compressible modes in the system. In this case the 
magnetic energy will cascade in parallel and perpendicular direction. One requirement 
is that magnetic and kinetic energy may cascade differently. The spectra suggest that 
this is indeed what happens in the simulations, but an explanation is hard. 

For (3 = 10 the two dimensional spectra are quite isotropic as we would predict it for 
a this widely hydrodynamical case (see figure [IT(cj| [TT(d)| [12(c)] [12(d)] [l3(c)| [l3(d)| 
I4(c71 and [11(d)] ). 



Similar investigations have been done by Vestuto et al. (2003) for the case of incom- 
pressible driving in the saturated state of the spectrum. He also found a perpendicular 
preference of the fluctuations for small values of the plasma-/?, which disappears for 
higher j5. 



5. Discussion 

In this paper we presented results from simulations of compressible MHD turbulence for 
low- and high-beta plasmas with compressible and incompressible driving. Our special 
focus was on the evolution of the turbulence. We found an obvious influence of the 
plasma-/3 and the type of driving on the anisotropy and the temporal energy evolution. 
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(a) /3=0.01 incompressible driving simu- (b) /3—0.01 compressible driving simu- 
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(c) /3=10 incompressible driving simu- (d) /3= 10 compressible driving simulated 
lated on a 512 3 grid on a 512 3 grid 



Figure 14: Pcomp in the saturated state of the turbulence simulated on a 256 3 and a 
512 3 grid. 
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Especially the influence of the driver has not yet been discussed for MHD plasmas so 
far. 

This may have influence on a number of astrophysical scenarios. When assuming su- 
pernovae as energy input, it is evident that e.g. the pile-up of fluctuation-energy due 
to slow tangling of magnetic field lines may be observed in the early stag lower 
temperature of the plasma. So by observing temperature evolution in turbulent driving 
it should be possible to determine which process is really driving the turbulence: Incom- 
pressible or compressible fluctuations. 

Even more striking is the influence on particle acceleration in evolving turbulence. Here 
also early supernova remnants are the object of interest. The first thing to observe here, 
is that a weak cascade in compressible driving may prevent low energetic particles to be 
accelerated due to missing energy at their respective resonant high wavenumbers. 
It is not completely clear if the anisotropy arising in first place will have an observable 
effect, but it is very well possible that the anisotropy in the magnetic fluctuations may 
change the transport and consecutively the acceleration and escape of cosmic ray from 
turbulent regions. 



6. Summary 

This paper gave summary of the evolution of turbulence in a driven MHD plasma. The 
evolution is described especially under the influence of different driving mechanisms 
(compressible and incompressible). First principles are able to explain the spectra for 
the very early development of turbulence. The anisotropy of wave generation is able 
to explain the specific features visible. For the further development of turbulence the 
plasma beta plays an essential role in the saturation of the spectra. We have shown the 
distribution of energy between magnetic and kinetic energy depending on the plasma 
beta. The results shown here are in agreement with earlier research done for hydrody- 
namics and goes into more detail regarding the specific anisotropy introduced with the 
background field. 
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A. Measure of the anisotropy 

Anistropy is one of the main features discussed in the context of this paper. The qual- 
itative behaviour of the anisotropy can be seen in the two-dimension plots of the spec- 
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Figure 15: Two-dimensional spectra for the sample spectrum Eq. 28 Shown are plots 
for the values of A = 0.1, 1, 10. 



trum. A quantitative description is far more difficult since the actual shape of the 
two-dimensional spectra may differ strongly. 

To give a rough quantitative comparison we are trying to boil down the plot to single 
numbers. We are using the contour lines to find a fe^ uso (A;|i). By following a contour line 
starting at a given value of fey from the ku axis to the k± axis we find this parameter. 

We want to illustrate the determination with a sample anisotropic function 



1 



(x 2 + Ay 2 ) 



12 



(28) 



The contour plots for different parameters A are shown in Fig. 15 For the simple case 
of A = 1 we find for example fc^ mso (2.3) = 2.3. The ratio of hi and k± is here 1, since 
the function is anisotropic. For A = 0.1 we find /c^ mso (1.7) = 2.7 and for A = 10 we find 

£80*0(2.6) = 1.6. 

This number gives a first quantitative hint on the anisotropy of the problem, but it is 
limited to axis aligned asymmetry. For feature propagating at 45 degree to the magnetic 
field it does not show any features. 
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